Convergence rate for the method of moments 
with linear closure relations 

Y. Bourgault 1 , D. Broizat 2 , P.-E. Jabin 3 

Abstract 

We study linear closure relations for the moments' method applied 
to simple kinetic equations. The equations are linear collisional models 
(velocity jump processes) which are well suited to this type of approx- 
imation. In this simplified, 1 dimensional setting, we are able to prove 
stability estimates for the method (with a kinetic interpretation by a 
BGK model). Moreover we are also able to obtain convergence rates 
which automatically increase with the smoothness of the initial data. 

1 Introduction 

1.1 Quick presentation of the moments' methods for 
kinetic equations 

We study here an unusual but very simple choice of closure for the moments' 
method where we can completely characterize the stability and convergence 
rates of the approximation. As far as we know this very simple situation was 
never considered before. Before presenting this choice though, let us very 
briefly give the main ideas behind the method of moments. 

Moments' methods have been introduced in [TO] in the context of the 
Boltzmann equation. This well known equation is posed on the density 
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f(t,x,v) of particles in the phase space and reads 

d t f + v-VJ = Q(f), (1.1) 

where Q is a non linear operator (in the velocity variable v) which expresses 
how the velocity of a particle may change when it has a random collision. 

Solving numerically an equation like fll.ip is in general very costly. The 
structure of the righthand side Q is non local in v and moreover the equation 
is posed in phase space which means that one has to work in dimension 2d+l 
if x, v G M. d (7 for instance if x G M 3 ). 

The moments' method is one possible answer to this problem and it con- 
sists in solving equations on the polynomial moments of /. Let m(v) be 
polynomial in v then 

d t <m,f > +V X - < v m(v), f >=< m, Q{f) >, (1.2) 

where we denote 

< m, f >= J m(v) f(t,x,v) dv. (1.3) 

Now instead of solving one equation in dimension 2d + 1, one has to solve 
several equations but in dimension d+1. Moreover in general one can expect 
that < m, Q(f) > is not too complicated to compute. 

However the system given by (11.21) is not closed as vm is always one degree 
higher than m. Therefore no matter how many moments mj one chooses, it 
is never possible to express all the v rrij in terms of the m^. This is the closure 
problem and it means that fll.ip is never equivalent to (11.2j) for any finite 
number of moments. 

Instead one typically chooses a closure equation, i.e. a relation between 
< vrrii, f > and the < rrij, f > for those i where vrrii cannot be expressed in 
terms of the rrij. 

The first big difficulty for this type of method is how to choose the clo- 
sure in order to ensure that the corresponding moments' system has good 
properties and gives a good approximation of Eq. fll.ip . This problem was of 
course recognized early on, see for instance [3], as well as the role of entropy, 
see [15] among many others. 

One of the first systematic ways of finding a closure was introduced in [12] 
and [13J. It is still not easy to actually compute the relation which means 
that it is often computed numerically instead (see [18], [6], [7]). Different 
closures can of course be used (see [2T] for example). 
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Theoretically even checking that the corresponding method leads to a 
hyperbolic system is not easy (we refer for instance to [I], [5]). Proving 
convergence rates seems to be out of reach for the time being although in 
practice it seems to be a good approximation (see [Hj for a numerical study). 

Let us also mention that the methods of moments has also been used for 
theoretical purposes (as in [S]) and not only numerical computations. 

We conclude this very brief overview by refering to [2J, [T7] or [T5] for 
more on numerical simulations for kinetic equations in this context. 



1.2 Linear closure relations 

The guiding question in this article is whether it can make sense to consider 
a linear closure relation. This is certainly delicate in the nonlinear case of 
Boltzmann eq. ( II. ip . Instead we choose a simplified Id setting where it is 
possible to fully analyze the method. 

Instead of fll.ip . we consider the linear model 

d t f + vd x f = L(f), (x,v)eRxR 

r , (1-4) 

Hf)= / Q(v,v*)f{t,x,v*)dv*-Xf 

with A > and where the operator Q corresponds to a velocity jump pro- 
cess. While much simplified with respect to f ll.ip . it is not uninteresting in 
itself, with applications to physics (see [8], [23]) or biology (see [H] for ex- 
ample). The equation had to be supplemented with some initial data, which 
for simplicity we assume to be compactly supported in velocity 

fit = 0,x,v) = f°(x,v) G L 2 (M 2 ), supp/°Clx/. (1.5) 

In general Q could even be assumed to depend on t and x. Here we make 
the additional approximation 




Q(v,v*) = q(v)y^ajV* ] l {( ^« )e/2} , (1.6) 



with q smooth and compactly supported in some interval I, d G N* and 
(otj)o<j<d ^ With this special form, one of course expects to be in a 
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very favorable situation for the method of moments. Hence this should be 
seen as a simple toy model where the method can easily be tested. 
Denoting the moments of the solution / by 

fi{(t,x) := ^v { f(t,x,v)dv, ieN, (1.7) 

Eq. (I1.4p simply becomes 

d 

dtf + vdj = L(f) = q(v) <*iA - X f- (1-8) 

3=0 

As we work in dimension 1, the structure of the hierarchy of equations on 
the moments is also very simple 

where we truncate at order N and we define the moments of q 

7i = n*= I v'qtydv, ieN. (1.10) 



In order to close the system, it would be necessary to be able to express 
in terms of the /if for i < N. The linear closure relation that we study here 
consists in assuming that is a linear combination of the lower moments. 
That means that instead of (\1.8\i or (|1.9|) . we solve 



dtHi + d x n i+ i = 7t I 2 ~ ^ h i = 0,...,N 



N 



1.11) 



i=0 



with (a , • • • , Oat) given real coefficients that one should choose in the "best" 
possible way. 

One can rewrite (II. lip in matrix form as 

d t M N + d x (AM N ) = BM N , (1.12) 
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where M N = M N (t,x) = {fi {t,x), . . . , fi N {t,x) f G R N+1 , and A,B e 
Mjv+i(M) are defined by 



/ 1 
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1.3 Basic properties of the method and the main result 

Given the simple form of (II. lip or (I1.12p . some properties are obvious. Given 
its linear nature, the system is hyperbolic if the characteristic polynomial 

N 

X a(X) = det(XI N+1 -A)= X N+l - (1.14) 

has N + 1 real roots, denoted by 

spec{A) = {X ,...,X N }. (1.15) 

This is enough to guarantee the well posedness of the numerical system (I1.12p 
but not necessarily good stability properties. Norms of the numerical appox- 
imation could for instance grows fast as iV increases. If the initial data is 
compactly supported in x then the solution is as well and the support prop- 
agates with speed max^ |Afc|. 

On the other hand, the major inconvenients of such a method are also 
pretty clear. For instance positivity of the even moments will likely not be 
preserved. Even the positivity of the macroscopic density fil has no reason 
to be propagated. 

However a careful analysis can in fact reveal that it is possible to choose 
appropriately the coefficients in order to have not only stability but also 
very fast convergence. 

For simplicity, assume that / = [—1, 1] (just rescale and translate other- 
wise) and choose the Tchebychev polynomial for \a or 

V0<fc<iV, \ k = cos (f^^\ n \ (1.16) 
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or 

*r^n(*--((f^>)). (i.i7) 

Then it is possible to show: 

Theorem 1.1 Assume that f° satisfies (|1.5p . i/iai Q satisfies (II. 6p . T/ien 

i/ie solution to the truncated moments ' hierarchy ( II. lip where the coefficients 
ai are chosen according to (I1.16P or (I1.17P satisfies 

sup sup ||^(R) <e TCd ' allqh2 ^\\f \\ L 2 {R 2 h (1.18) 

0<i<Artg[o,T] 

1 /2 

where Cd, a = y/Md + 1) (j^Lo ^ 2 ) ^ s independent of N . 

In addition if f° G H k (M. 2 ), d = 0, denoting by f the corresponding 
solution to (11.41) and defining its moments by (11.71) . one has 

C 

ll^o - ^olU°°([o,r],L2(]R)) < Nk _ 3/4 x ||/° || H fc( K 2). (1.19) 
where C > depends on T,X,q and k. 

Remark 1.2 The convergence result is given for the first moment and d = 0. 
Higher moments would lose small powers of N so the same proof would give 

/II CN k{i+d)/N 
\W - Hi ||£°°([0,T],£2(R)) < Jjk^i X \\j ll^ fc (K 2 )' 

We do not know whether those results are optimal or not and in particular 
whether the numerical value 3/4 in (11.191) is optimal (though actual numerical 
simulations do suggest it is). 

Hence in this setting, the conclusion of this analysis is that the linear method 
of moments should be seen in the same light as spectral methods (see [16] for 
instance). It is stable, converges automatically at order k — 3/4 if the initial 
data is in H k but does not propagate any additional property (positivity 
being probably the most important). 

The next section is a more detailed presentation of the stability and 
convergence analysis in the general case (without necessarily choosing the 
Tchebychev points). The corresponding technical justifications and calcula- 
tions are presented in two separate sections. Theorem [1J] is proved in section 
5. The last section is an appendix and recalls some well known results. 
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2 Stability and convergence results 



We present here in more details the kind of stability and convergence results 
that can be proved for the method (11. lip without assuming any particular 
choice of the eigenvalues (like (11. 16ft ). We give here the main ideas in the 
approach and leave the technical proofs to further sections. 



2.1 Eigenvectors for System (11.121 ) 



As they enter in some of the estimates, we start by a short parenthesis about 
the eigenvectors for the matrix A. Define P G Mjv+i(K.) the matrix of the 
eigenvectors; then 

Pi. j A}. 0<z,j<iV. (2.1) 
Its inverse P~ l can be computed as easily and 

Prl = —, 0<i,j<N, (2.2) 



7T, 



where 



7T< 



H^-A,), /<„ ( n v j £ W^^H^rm- ( 2 - 3 ) 

jyti k 1 <...<k N _ j 1 = 1 l = Q \ 

with the convention that pi^ = 1. 

Moreover, an easy computation shows that: 

= K ~ J - °ArAf a j+2 \i-a j+1 , < z, j < N. (2.4) 

We can notice that pij is an homogeneous polynomial of degree N — j in the 
eigenvalues (A , . . . , Ajv). 

2.2 Stability estimate and kinetic interpretation of the 
method 

Let us first state our main stability estimate. 
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Theorem 2.1 Assume ( II. 5p and that {X , . . . , A^} C J. Moreover assume 
that there exists a function Pn{v), positive on I, such that 

f ^p-dv = V>(A fc ), Vi? G M 2W+ i[X]. (2.5) 

TTien, t/ie hyperbolic system (11.12p stable and 

/ N \ V2 
\fc=0 / 

(2.6) 



sup ||//i(t)||jy»(R) < 
te[Q,T] 



where 



C N (f°) =( ff \f°(x,v)\ 2 p N (v)dxdv) ' , (2.7) 
C N (q) = A N4 ( [ \q(v)\ 2 p N (v)dv) ' - A, (2.8) 



A 



N,d 



N 



j=0 k=0 



N 

£ 

fc=0 



^2d+2 



(2.9) 



This result does not assume any particular distribution on the eigenvalues 
but of course it is by no means guaranteed in general that one could find p^ 
satisfying ( 12. 5p . Notice that the corresponding relation is really a quadrature 
formula for computing integrals on / which we ask to be exact for polynomials 
of degree 2N + 1. 

We do not prove this result directly on the system (I1.12p . Instead we show 
a corresponding result on a linear BGK problem (see [1] for the simplification 
of collision kernels into BGK models). 

For any e > and N > d, consider the equation 



d t f £ N + vd x l 
fff(0,x,v) = f (x,v 



(2.10) 
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where the linear operator L is defined by (jl.8p and the Maxwellian M^ N > : 
f i — y M^ N 'f satisfies the moment conditions 

v { M {N) fdv = [ v l fdv i = 0,...,N, 

i J i 



f v N+i M Wf dv = V fl! [ v i fd . 
Ji i=0 Ji 



This problem is a kinetic approximation of the macroscopic problem ( II. lip 
as formally when e — > then one recovers (jl.llj) from (12.1 Op . 

Indeed, multiplying (I2.10p by v % and integrating over v G /, we obtain 

9tl4 + 9 x f4 +1 = ji - for * = 0, ... , N, 

where \i\ := J I v l ff f dv for % 6 N. Moreover, when e — > 0, we have formally 

N r N 

e _ / „.JV+i re j„, _ . / „JV+i »>f(JV) fe i„, _ „ / „,i J?e j„, _ \ r „ ,,£ 



which is our closure relation. 

The interest of (I2.10p is to make some computations more transparent 
and easier to follow and in addition at the kinetic level which is more natural 
given the original equation. If we can prove stability estimates for (I2.10p that 
are uniform in e then simply by passing to the limit, we will obtain estimates 

for fln2j-flnu. 

The most obvious choice for the maxwellian is simply 

M {N) f = ]T ( I v'fdv) m u (2.12) 
i=o ^ I ' 

where the v i— > rrii(v), i = 0, . . . , N, are any functions satisfying 

v j mi{v)dv = 5ij, 0<j<N, ^v N+1 mi(v)dv = a { . (2.13) 

The conditions (12. 13[) ensure that (12. lip holds. 
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There are obviously many ways to choose the raj s.t. (I2.13P is satisfied. 
What we are looking for is a choice compatible with an inner product $tv 
such that the application is an orthogonal projection for $at. Formally 

this implies that for any / 



In addition this inner product should have good symmetry properties s.t. for 
/ and g 

$N(f,vg) = $ N (vf,g). 
The simplest way to ensure this is to look for a weight p^ s.t. 



This is the strategy that we implement. Find appropriate conditions on p^ 
and the rrii to obtain the correct structure and then simply bound $jv(/, L(f)) 
in terms of <!>#(/)■ 

For the first part of this strategy we actually prove 

Theorem 2.2 Let Pn{v) a positive function on I such that 



$ N (f,MWf)<$ N (f,f). 




In that case formally 



$>N{f,vd x f) = -d x <t> N {f,vf). 



Therefore if fff solves (12.101) then one expects that 






(2.14) 



We set 




Then, the map 4>n defined by 





(2.15) 
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is an inner product on E^, and the Maxwellian M^ N > : — > E^, defined 
by 

is an orthogonal projection and satifies (12. lip , where (Tj)o<i<Ar is the basis 
ofM.N[X] defined by: 

N 

f t (X) = Y,Qk tl X k , 0<t<N. 

k=0 

with Q = {P^P- 1 . 

Furthermore, with that choice of M^ N \ any solution f N of the problem 
(12.1 Op formally satisfies: 

^ / Mf £ N,fN)dx<2 [ MfN,L(f £ N ))dx. (2.16) 
at Jr Jr 



2.3 Error estimate 

We now turn to the convergence of the method of moments to the solution. 
For simplicity we state the results here for the case d = in (II. 6p . namely 

Q (v,v*) = 

We also assume that / C [—1, 1], still for simplicity. 

The convergence results of course require a stability estimate. However 
in themselves, they do not use the specific form of the closure. 

Therefore here we do not assume any specific closure relation. Instead 
we assume that we have a well defined methods of moments, i.e., some way 
of computing /ij which satisfy 



dtHi + d x fi i+l = 7i/i - Xfii, i = 0,...,N 
Hn+i = F{ni, [in). 

where 

li = / v l q(v)dv. 



(2.17) 
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Moreover we assume that the corresponding method has good stability esti- 
mates in the sense that 

3C, 7 > 0, For any (fii) i= o.. N and (fii)i= ..N solutions to (11. lip for /°, /°, 

N 

then ^2 WfM ~ frill l°°{[o,t], H k (R)) < CN^Wf — / ||l2(/, H k (R)), 

8=0 

(2.18) 

where C, 7 may depend on T, q, A but not on N or f°. 

Note that Th. 12.11 can indeed be expected to imply such a result. The 
exponent 7 will depend on the choice of the coefficients in our method. Of 
course Th. 12. II controls only the L 2 norm of the /ij. However as the method 
fll. lip is linear, the are also a solution to the same system and an 
estimate like (I2.18P can be derived. For a more detailed analysis of how to 
obtain (12.18P from Th. I2.1[ we refer to Section [5] where it is performed when 
the Afc are the Tchebychev points. 

For any method that satisfies (I2.18p . then one has the following conver- 
gence result 

Theorem 2.3 Assume ( II. 5p with I C [— 1, 1], that the method (12.17P sat- 
isfies (I2.18P for some 7 > and that f° G H k (R 2 ), with fcef. Consider f 
the solution to ( II. 4p with d = and the corresponding solution ^ to ( I2.17P . 
Then, for all T > and for all N > 1, we have the estimate 

C 

11^0 - ^olU°°([0,T],L 2 (R)) < ^ fc _ 7 X ||/°||fl-*(K a )- (2.19) 

where C > depends on T,X,q and k but not on N. 

This error estimate is a sort of interpolation between the stability bounds 
and the following result for C°° solution to (12.171) 

Proposition 2.4 Assume that the (fJ<i)i=o...N+i solve 

dtut + d x /j i+1 = 7^0 - A//i, i = 0, JV, (2.20) 

with fii{t = 0) = for i = 0...N. Then, for all T > and for all N > 1, we 
have the estimate 

rpN N 

II M0 - MolU°°([0,T],L2(R)) < C X — ll^ r Mi|U°°([0,r] 1 L 2 (R)), (2.21) 

i=0 

where C > depends on T, A, q. 



12 



3 Proof of Theorems 12.11 and 12.2 



To simplify the presentation, we omit here the subscript iV in En, 4>n, Pn, 
and the superscript N in M^ N \ 

3.1 Elementary space decomposition 

The difficulty is to combine the fact that M has to be an orthogonal projec- 
tion for $ with the symmetry property on $. We take here a slighty more 
general approach by not assuming directly that $ satisfies (12.151) . 

Consider a Maxwellian which has the form (I2.12p and satisfies (|2.13|) . 
First, such an application M is a projection because M o M = M, which is 
a straightforward consequence of 

M(rrij) = I / v^rriidv J rrij = rrii, i = 0, . . . , N. 

Moreover, one has that 



KerM=^feE, J v l f(v)dv = 0, i = 0, . . . , Nj := K, 

Ker(M — I) = Im(M) = Span(m , m 1? . . . , tun) '■= V, 
and we have the space decomposition 

E = V © K, 

with 

dim(V) = codim(K) = N + 1. 

We start by a more detailed explanation of the sufficient conditions to obtain 
Th. 12721 

Lemma 3.1 Assume that the inner product : E x E — >■ K. satisfies 
K = V (i.e. the decomposition V © K becomes orthogonal). 

V(f,g)eE\ <Kvf,g) = <l>{vg,f). (3.1) 
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Then, for any solution = f^(t,x,v) of the problem (12.101) . inequality 
(I2.16P formally holds: 

/ ||^(t,x,.)||}dx<2 / <P(f N ,L(f N ))dx. 
ai Jr. Jr 

Proof : 

Take a smooth / = f(t, x, v) such that d t f + vdj = \{Mf - f) + L(f). 
We have 



/ \\f(t,x,.)\\ldx = 2 [ <j>{f,d t f)dx 
at Jr Jr 



- [ 4>(f,Mf-f)dx- [ <j>(f, vd x f)dx+ [ 4>(f,L(f))dx 
e Jr Jr Jr 

Since Mf - f 6 K and Mf e V, we have (f>(Mf, Mf - f) = 0, thus 

<P(f, Mf — f) — <j>(f - Mf, Mf-f) = -||/ - M/|U < 0. 
We deduce 

^ / \\f(t,x)\\ldx<-2 [ <f>(f,vd x f)dx + 2 [ <P(f,L(f))dx. 
at Jr Jr Jr 

Thus, having J R (f>(f,vd x f)dx = is sufficient to obtain ( 12. 16ft . 
Since <f>(f,vd x f) = d x (j)(f,vf) - <p(d x f,vf), we can write 

4>(f, vd x f)dx = - I (4>(f, vd x f) - <j>(vf, d x f)) dx. 
1 Jr 

Therefore the symmetry condition (13.11) on <fi is enough to conclude. If / is 
not smooth enough to follow the previous steps, one simply regularizes it by 
convolution in i. As the equation is linear and x is only a parameter in M 
and L then the regularized function solves the same equation. Therefore it 
satisfies ( 12.161) and letting the regularizing parameter vanish, one recovers 
the same inequality for /. □ 

Now, take a inner product <fi on E such that K = V^^ , i.e. K and V are 
orthogonal for this inner product. The symmetry condition (13. ip is obviously 
equivalent to 

V(f,g)EV 2 , <P(vf,g) = ^vg,f), (3.2) 

V(/, g) e K x V, cf>(vf, g) = <f>(vg, f), (3.3) 

V(f,g)eK 2 , <Kvf,g) = <f>(vg,f). (3.4) 
We study each of those conditions in the following subsections. 
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3.2 Study of the condition flgg ) 
Proposition 3.2 T/ie condition (13.21) zs equivalent to 

A T Q = QA, 

where A is the matrix defined by HI . 1 3[) and Q is the symmetric definite 
positive matrix defined by 

Qij = <f>(mi,m,j), 0<i,j<N. 

Proof : The condition (I3.2p is equivalent to 

(p(vrrii, rrij) = (j)(rrii,vmj), VO < i,j < N. 

Let G {0, . . . , N} 2 . Since rrij G V, he have 

4>(vmi,mj) = (f)(M (villi), rrij). 

Moreover, (12.131) implies 

M(vrrii) = nii_i + aim N , < i < N, (3.5) 

with the convention m_i = 0. Thus 

<p(vm u rrij) = <f>{mi-i, m j) + a>i4>( m N, rrij) = Qi-i,j + <hQNj- 

But 

N 

(A T Q)ij = Ak t iQk,j = Qi-i,j + a-iQNj, 

with the convention = 0. Therefore, we have 

(f)(vmi,rrij) = (A T Q) itj . 

Thus, (13. 2p amounts to the matrix A T Q be (real and) symmetric, i.e. A T Q = 
QA, since Q is symmetric. □ 

This result suggests a particular way of defining the inner product on V: 
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Corollary 3.3 Choose s.t. 

Q = {p -l ) Tp-l = { p p T ) -l ; 

where P is defined by (12. ip . This choice implies 



N 



(m l ,m J ) = Q h3 = J2PkjPkj- (3-6) 



fc=0 



Proof. We have, denoting by D = diag(Xo, . . . , Ajv): 

A T Q =A T (P- 1 ) T P- 1 = (P- 1 A) T P- 1 = (DP-yp- 1 = (p-YiDP- 1 ) 

= (p-yp^A = QA, 

thus Q satisfies A T Q = QA. Moreover, it is easy to check that Q is symmet- 
ric, definite, and positive. □ 

In the rest of the proof, we always choose $ according to (13.61) . 



3.3 Study of the condition (13.31 ) 



Now, we assume (13. 6 1) to be satisfied and analyze ( 13.3H . 

Proposition 3.4 Assume f 1 3 . 6 1) . and consider the following polynomials 

1 N 

Ti{X) = 7 ^Y J Qk,X\ < i < N, (3.7) 

Qo,n k=Q 

where the matrix Q G M^v+i (R) is defined by (13. 3|) . 
Then, the condition (13.31) is equivalent to 

T N (v)mi(v) = Ti{v)m N (v), v el, < % < N - 1, (3.8) 

w , _ v , ( 1 v XA(v)m7v(f) ,\ f N+lf/ v, , v 

V/ G if, - — x — — — , f \ = v + f(v)dv. 3.9) 

Wo, at -Tjv(«) / 7/ 

The proof of this proposition is split in several lemmas. First we find two 
equivalent conditions to (13.31) . 
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Lemma 3.5 The condition ( 13. 3 j) equivalent to 

Vi G {1, . . . , iV}, t>mj — mj_i — ajmAr = ^ ,JV (vm — a m N ) (3.10) 

Qo,n 

V/ G K, <f> ]K (-^—(vmo - a m N ), f) = [ w N+1 f(w)dw, (3.11) 

\WO,N J J I 

We then have to study conditions (I3.10p and (13. lip . 

Proof of Lemma 13.51 The condition ( 13. 3 j) is obviously equivalent to 

<f>{vm i ,f) = <f>{m i ,vf), 0<t<N, f e K. 

Let i G {0, . . . , N} and / G K. Since f & K, he have 

0(wii, /) = (j)\ K {vmi - M(vmi), f) = <\>\K(vmi - m;_i - a^m^, /), 

using (I3.5p . Similarly, since mj G V, 



(m t ,vf) = <j>\ v {m h M{vf)) = Q i>N I w N+1 f{w)d 



7 



where the last equality is deduced from 

V/ G K, M(vf) = ^w N+1 f(w)dwj m N , 

recalling that j w n f dv = 0, Vn < N. 

Thus, the condition (I3.3P amounts to for any < i < N and / G K 

(j)\ K (vmi - rrii-i - aiin N , f) = Q i)N J w N+1 f(w)dw. 

from which we deduce (13. lip , and 

\7^—( vm i - m i-i ~~ a i m N) ~ -^—(vm - a m N ), A = 0. 

Therefore, for alH G {1, . . . , N}, we have 

I 1 

— — (vrrii - m^i - aim N ) - — — (vm - a m N ) G K n = {0}, 

which shows the relation (I3.10p . 

Conversely, the conditions ( I3.10p + (l3.1ip imply ( 13 .3p just by following the 
previous steps in reverse order. □ 
We start with condition ( I3.10p 
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Lemma 3.6 Consider the polynomial 



N 



8=0 

where 

Qo,N <P\jriQ, m N ) 
The condition (13.101) is equivalent to 

D(v) mi (v) = Ti(v)m N {v), v el, 0<i<N, 

where (Tj) <i<7v are the following polynomials: 

T N (X) = D{X), 

T^X) = XT t {X) - ai D(X) - iX A(X), l<t<N. 



(3.12) 



Proof. Setting ft = for i e {0, . . . , N}, we deduce from fl3TTU|) the 
recursive formula 

m ._ 1 = vrm - (at - (3ia )m N - ft(vm ), i = 1 . . . , N, 

which leads to, for any < i < N 

tiV-i-l \ /N-i 
v N-i _ v k (a i+k+ i - /3 i+k+ ia ) J m N - I ^ f3 i+k v k | /»o- 

fc=0 / \k=l 

Thus 

AT-l 



D(w)m (f) = ( v N - ^2(a k+1 - f3 k+1 a )v k ] m N (v), 

where 



fc=0 

D(v) :=p + p lV + ---+p N v N 

We deduce 



vD(v)m Q (v) = (xa(v) + a D(v))m N (v). 
Thus, it comes from the recursive formula on the m ; that 

D(v)rrii-i(v) = vD(v)rrii(v) - {aiD(v) + (3iXA(v))m N (v), i = 1 . . . , N, 
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which is exactly 

D(v)m i (v) = T i (v)m f f{v), i = 0, . . . , N, (3.13) 

with 

T N (X)=D(X), 

T^X) = XT t (X) - a t D(X) - /3 iX a(X), i = l,...N. 

Conversely, if the functions (?7ij)o<i<jv satisfy these relations, then (13.101) is 
obvious, almost everywhere v G / (except at the roots of the polynomial 
D). n 

Remark 3.7 We notice that 

XT (X) - a D(X) = X a(X). (3.14) 
Furthermore, the formula (13.121) implies: 

1 + deg(Ti) < max{deg(T i _ 1 ), N + 1), 1 < i < N. 
But T G R N [X], thus Ti G R N [X] for all i & {0, ... , N}. 
We can now prove 

Lemma 3.8 We have the explicit formula (13.71) for the (Tj) <i<v which is 
recalled here 

1 N 

Ti(X) = —— V Q k>i X k , 0<i<N. 
Qo > N k=0 

Proof : First, according to the definition of D in Lemma \3. 6 1 we have, for all 

ke{0,...,N}, 

^ ^ V V 

1 

Qo,N T^k 
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Then, the recursive formula (13.121) easily implies: 

Tj{\k) = (A^ ~ J - oatA^^ 1 a j+2 \ k - a j+1 )D(X k ), < j < N 

which gives (according to (I2.2p - fl2.4l) ) 

Tj(\ k ) = Pk,jD{\ k ) = J k ' 3 = < j, k < N. (3.15) 

Since each polynomial T{ have a degree equal to N, it is enough to check 
equality (13. 7\i on the set {Ao, . . . , Xn}, which allows to conclude. □ 

Remark 3.9 The formula (13. 7p shows that the polynomials (Tj)o<i<Ar are a 
basis o/M.n[X] since the matrix Q is invertible. 

We may finally characterize (13. lip 

Lemma 3.10 We assume that (I3.10p holds. Then, the condition (13.111) is 
equivalent to 

uf r- v x ( 1 XA(v)m N {v) \ f N+1 

Vf EK, <j)\ K [ - — x — — — , / = / v N+l f (v )dv. 

\Qo,N T N (v) J J I 

Proof : It is straightforward using (13. 13[) and the formula (I3.14p . □ 

We now have all what is needed to prove Prop. 13.41 
Proof of Prop. 13.41 By Lemma I3.5[ condition (13 .3p is equivalent to (I3.10P - 
(13. lip . By Lemma IBTBl condition (13.101) is equivalent to (13 .8p provided that 
the Tj are defined by (13.121) . Lemma 13.81 shows that the recursive formula 
(I3.12p actually gives the explicit formula (13. 7p . Finally by Lemma I3.10[ we 
know that condition (13. lip is equivalent to (13. 9p thus concluding the proof. □ 



3.4 Study of the condition (J23D 

We prove 

Proposition 3.11 Assume (JSSD+fJXHD+fJSS])- Then, setting 

, x Qo,nTn(v) 

P[v) := pr— , (3.16) 

m N {v) 

and 

Mf, 9) = J i f(v)g(v)p(v)dv, V(/, g) e K 2 , 
the condition (13. 4p is satisfied. 
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Proof : This choice is compatible with (13. lip , since (13. lip can also be written 
as 

V/ € K, <j> [K {^j, /) = Jv N+1 f{v)dv, 

and since we have JjXA^f^dv = JjV N+1 f(v)dv since / e K. 
Moreover, (13. 4p is obviously satisfied with this choice. □ 



3.5 Choice of on the subspace V 

For the moment $ is defined as a weighted L 2 type inner product on K by 
Prop. 13. Ill and on V by Corollary 13.31 

We wish to define $ as a weighted inner product on the whole K © V . 
The following lemma shows that provided p satisfies the right relations then 
$ as defined by Prop. 13.111 and Corollary 13.31 is automatically of the right 
form. 

Lemma 3.12 Assume (I3.6p . (13. 8p . ( 13 .9p . (I3.16p . and assume that the weight 
p satisfies the moment conditions: 




dv = J2 R ( X k), R€R 2N [X]. 



Then, we have: 

<l>\v(f>9) = Jf(v)g(v)p{v)dv. (3.17) 



Proof : It is sufficient to show that the formula ( 13. 17ft holds for (/, g) 

dv 



(mi, rrij). We have, for e {0, . . . , N} 2 , 



mi(v)mj(v)p(v)dv = / (fifjj (v) x 



p(v) 



setting 



ft = Q , N T i} < % < N. 

Moreover, we have 

N N 

<f>( mi , mj ) = Qu = j2 p k~J p kj = E (^) ( Afc )' 

k=0 k=0 
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according to (13. 6p and (I3.15p . 

Since the ( JV+1 K Af + 2 ) polynomials {TiTj)Q<i<j< N are in R 2 ivM> we see that 
the assumption on p guarantees that (13.171) is satisfied. □ 

3.6 Proof of the theorem (12.21 ): Synthesis 

We summarize here all the definitions and check rigorously that they are 
compatible. 

So, assume there exists a function p = p(v), positive on / and satisfying 
the moment conditions: 

f ?^l dv = J2R(\ k ), Vi? G R 2 n+i[X]. (3.18) 

Note that this in particular implies that 1/p is integrable on I. 

The space E = L 2 (I, p{v)dv) is a Hilbert space, for the inner product 



<Kf, 9) = J f(v)g(v)p(v)dv, (f, g) G E 2 . 
We define the map M : E — » E by 



TAv 



i=0 ^ I 



p{v) 



where 

N 



f l {X) = Y,Qk^X\ 0<i<N, 



fc=0 



and Q = (P T )- 1 P^ 1 (the matrix P is defined by (l2TTj) ). 

• First, the map M is well defined as V/ G E, Vz G {0, . . . , iV} 



v % f(v)dv 



I |2i \ 1/2 / » ^ 1/2 

< ( / / l/(«)lV(«)cfo 



and thus M/ makes sense. Moreover, for / G E, we have M/ G 
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because 



JjMf(v)\ 2 p(v)dv = J 



N 



E ( f/f^) ^ 



i=0 



p(v)dv 



< 



N 

E 



,i=0 



N 

<'E 



w l f(w)dw 



w l f(w)dw 



N 

E 



i=0 



X 



N 

E 

i=0 



p{v)dv 



dv < oo. 



• It is obvious that M is linear. We check that M is a projection: let 
f e E,we have, for < j < TV, 



J^v j Mf(v)dv = (J 

w l f(w)dw J / 



i=0 
N 



v j f t {v) 
p(v) 



dv 



i=0 ^ 7 ' \fe=0 / 



Moreover, we have 

\{ = P jtk , T,(A fe ) = Prl, < i,j, k < N, 



thus 



We deduce 



A? 



52\if i (\ k ) = 8 iJ , 0<i,j<N. 



k=0 



J^v j Mf(v)dv = Jw j f{w)dw ) 0<j<N, 
which implies M o M = M. 

M is an orthogonal projection (for the inner product 0), because it is 
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a self-adjoint projector: 
<f>(Mf,g) = [ Mf{v)g{v)p{v)dv 



1 

N 



= E (^7H^) (^ff t {v)g{v)dv 

= EE^ ( / <fH<^) ( / A(«)*>) = <P(f,Mg), 

as the matrix Q is symmetric. 

The Maxwellian M satisfies moment conditions: 

tfMftydv = ftff^dv i = 0,...,N, 

i J i 

[v N+1 Mf(v)dv = J2ai ftffWdv 
Jj i=0 Ji 

In fact, the first conditions have been already established, and the 
second ones result from the following formula, using (13. 18[) and since 
XA{v)Ti{v) is a polynomial of degree 2N + 1, 



i=0 

N , „ \ TV 



( f wj( w )dw) j2xA(x k )m k ) = o. 



• By Corollary 13.31 Prop. 13.111 and Lemma 13.121 the inner product $ 
satisfies the assumptions of Lemma 13.11 This concludes the proof of 
Theorem 12.21 □ 

3.7 From Th. E21 to Th. EUl Uniform stability esti- 
mate on the BGK model 

The first point is to control the collision term which is done by 
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Lemma 3.13 Define p and $ as per Th. \2.2\ then 



j^jj £ N L{f N )p{v)dvdx < (^A N)d (^j\q(v)\ 2 p{v)dv 

\fx\ 2 p(v)dvdx, 



1/2 



A 




(3.19) 



with 



A 



NA 



\ i=0 \ 3=0 k=0 \ j=0 \ k=0 



d N 



EE*? 



E 



N 



1 \ 2d+2 

. 1-A 2 ' 



Proof: We have 




fx(t, x, v)L(fx)(t, x, v)p(v)dvdx 



So 




/ / fff(t, x, v)q(v) V f / ^'/, 
Jr J i =0 W/ 



^ / / l/^(M, u )l 2 p( u )dwfc- 



/Ir(*> x ; v ) L (fN)(t, x , v)p(v)dvdx 



'ff(t,x,w)dw I p(v)dvdx 




\ j= o Jl 



w 3 f e N {t,x,w)dw\ \ f e N (t,x,v)q(v)p(v)dv 



- A / / \fx(t,x, v)\ 2 p(v)dvdx. 
Jr J i 

We can control the moments of in the following way: 



w J ffj(t, x, w)dw 



and the assumption f)2.14p implies 



w 3 f^(t, x, w)dw 



N 



,fc=0 
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We deduce that 
d 



j=0 



Moreover, we have 



d N 



£ ED? w.a>' 



,j=0 £i=0 



f £ N {t,x,v)q(v)p(v)dv 



<P(fN,qy<MQ,Q)<P(fNJ e N)- 



Thus, combining the last two inequalities and using Cauchy-Schwarz inequal- 
ity, we obtain 



f £ N L (f £ N)p( v ) dvdx <* 



\ 3=0 \ 3=0 



d N 



k=0 

-A / (p(f £ N ,f £ N )dx, 
Jr 

which is the desired estimate. □ 

Combining Th. 12.21 and Lemma I3.13[ we can obtain stability estimates 
for (12.101) uniform in e: 



Proposition 3.14 Let f e N a solution of (I2.10p with f ll .5 j) . Then, the follow- 
ing estimate holds for any t > 

sup/ \f £ N {t,x,v)\ 2 p{v)dvdx <e tCN ^ [ \f°(x,v)\ 2 p(v)dvdx, (3.20) 

e>0 Jr J i Jr J i 



where 



and we recall 



1/2 

Cv.. / = 2(A.v., / ( j \q(v)\ 2 p(v)dv\ — A 



A 



N,d 



d N 

EE* 

j=0 k=0 



N 

E 

k=0 



1 - \f +2 

1-A 2 , ' 
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Proof : It is straightforward using Gronwall lemma, as according to (I2.16P 
and (I3.19p . we have 

\fx(t,x, v)\ 2 p(v)dvdx. 




i 



□ 

We are now ready to conclude the proof of Th. 12.11 Now, we show that 
we can pass to the limit e — > in the BGK model (12.10)) to obtain a stability 
estimate on the hyperbolic system (I1.12p . 

We fix N > d. Using (13.20)) . we can see that the family (fff) e >o is bounded 
in the space L 2 (]0, +oo[i oc xM, x xI v , p{v)dtdxdv). Thus there exists a sequence 
Ek — > such that 

fc— ^oo 

f £ k v f., r- T 2 



iocXlRj. x I V: pN^v^dtdxdv^ . 

k— >oo 



Therefore 



d t f N k +vd x f^ - d t f N +vd x f N , L(f N k ) - L(f N ) in VQO, oo[xlx/), 

k— >oo k— >oo 

from which we deduce 

MUn) ~ fN = £k(d t f N h + vd x f N k - L(f N k )) - in V'QO, oo[xR x /), 
and thus 

M(f N k ) - f N in V'()0, oo[xR x I). 

Hence, passing to the limit in (12.111) . we show that the function is a kinetic 
interpretation of the system (j!.12p . in the sense that (/ij = J I vl fNdv) Q< . <N 
satisfies (|1.12p . As this system is linear, hyperbolic, it has a unique solu- 
tion for a given initial data. That means that any solution of (I1.12p can 
consequently be obtained as the moments of a limit of ffj. 
Moreover, fx also satisfies the bound (13.201) : 

Vt>0, f [ \f N (t,x,v)\ 2 p(v)dvdx < e tCN >i [ [ \f°(x,v)\ 2 p(v)dvdx. 

JR J I JR J I 
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Thus we obtain the estimate (12. 6p . since 



\lH(t,x)\' 



v l fN(t, x, v)dv 



< (l>?) f\Mt, 

\k=0 J Jl 



x, v)\ 2 p(v)dv, 



according to the assumptions on p^. The proof of Th. 12. 11 is now complete. □ 



4 Error estimate: Proof of Th. [273] and Prop. 



2.4 



4.1 Estimates provided by the model 

In order to prove Th. 12.31 we will need good smoothness properties on the 
solution to the exact equation (11.4j) . Fortunately this model is very simple 
to manipulate and the estimates we need easy to obtain in that case. 
Let us first start with the support in velocity 



Proposition 4.1 Assume (II. 5ft . (II. 6ft with I C [—1, 1]. Then the solution 
to fll.4p satisfies 



Vt > 0, 



a.e. x G 



supp v f(t, x, .) C I C [-1, 1]. 



Proof : If / is a solution of (ll.4p with d — 0, then, using the Stokes formula, 
we have, formally 

^ J \f{t 1 x ) v)\ 2 dx = 2 y f(t,x,v)q(v) (^j f{t,x,v*)dv*^jdx 



-2A / |/(t,x,v)| 2 cfe 
Jr 



<2g( W ) / /(t,a;,t;)( ^f(t,x,v*)dv*)dx. 



Therefore, integrating in time, we obtain 
\f(t,x,v)\ 2 dx< [ \f°(x,v)\ 2 dx 



+ 2q(v) / / f(s,x, v) I / f(s,x, v*)dv* ) <£eg?s, 

and since suppg C /, we get the result. □ 

The model (ll.4p also propagates the if fc -norm of the solution 
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Proposition 4.2 Assume (ll.5p . (II. 6p with I C [—1, 1]. Then the solution 
to (II. 4p satisfies 

vt>o, ||/(t)|| 

where C = \\q\\L 2 (i) — A. 

Proof : First note that for any k, d^f is also a solution to (II. 4p with d%f° as 
initial data. Then 



^/(t,x,t;)g(t;)^ l 8PJ(t,x,v*)dv 



-2A||/(*)|| 

Using Prop. (14.11) and Holder inequality, we easily obtain 

^ll/(*)llfr*(R, l 2 (/)) ^ 2 (IMU 2 0O - A ) ll/(*)llfr*(H, l 2 (/))> 

and a simple Gronwall lemma gives the result. □ 
We may finally conclude from Props. I4.1H4.2I that 

Corollary 4.3 Assume (11.51) . (11.61) with I C [—1, 1]. Then the moments of 
solution to (11.41) satisfy 

V*GN, Vt>0, ||^ / (*)|U 2(R) <e Ct ||/°|| z2(K2) , (4.1) 

and 

N 

Eikwiu^ (2iv+i) i /2 e «|| / o| U2(M2) . 

where C = \\q\\L 2 (i) — A. 
Proof. Notice that 

|/if I < J Jv\ i f(t,x,v)dv< ^_L_ (J Jf(t,x,v)\ 2 dv^ 
by Cauchy-Schwarz. Hence 



|^(t)|U 2 (M)<^==||/(t)||L 2 (M 2 ), 

+ t 



and one concludes by using Prop. 14.21 



□ 
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4.2 Proof of Prop. 12.41 : Error estimate in the smooth 
case 

Since the functions are smooth in the space variable, we have by (|2.17j) . 
for all < i < N, 

d f 

-77 llAiiC*) IliariR) = 2 / Hi(t,x)d t Hi(t,x)dx 

'<' JR 

= -2 / [ii(t,x)d x iJ,i +1 (t,x)dx , 



+ 2% / x)/io(t, x)dx — 2A / \[ii(t, x)\ 2 dx, 
Jr Jr 

or 

^H/^(*)lli2( R ) < 2 II Vi(t) II L 2 (M) II d x p i+1 {t) || L 2 (R ) 

+ 2|7<| ||^(t)|| L 2 (R) ||/i (t)|| L 2 (R) - 2A||^(t)||| 3(R) . 

Thus, we obtain 

^11/^0) Hl 2 (ir) < ||^i+i(*)|U2(R) + |7t| IM*)||z,2(r) - A||/ii(t)|| L 2 (R) . 

Of course the same computation can be performed on the d k [ii, obtaining 
VA;GN, Vie{0,...,N}, 

l\\d k xf , t (t)\\L H R) < \\d k x +1 » i+1 (t)\\ Lm + N RVo(*)|b W (4-2) 

- \\\d k x fIi(t)\\ L 2 m . 

This sequence of inequalities lets us control ||/i (t)||,L 2 (]R) in term of the "last 
derivatives", namely (||^A t i||L 2 (R)) <i<j\r To do that, we set 

k 

Vt>0, Vke{0,...,N}, H k (t):=J2\\d k x ^h'HR)- 

i=0 

The coefficients 7« are easily bounded by H^Hl 2 - Hence the inequalities (14.21) 
imply 

Vt > 0, Vfc G {0, . . . , N}, j t H k (t) < H k+1 (t) + (2||g|| L 2 (/) - A) H k (t). 

(4.3) 

We now use 
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Lemma 4.4 Let (H k (t))k>i be a sequence of nonnegative C 1 functions such 
that: 

v* eN , vt>o, (flfif>(')+M() , 

where C > is a numerical constant independant of k. 
Then, we have: 

Vp G FT, Vt > 0, H (t) < fit - sf-^c^H^ds. 

Proof : First, the assumption may be rewritten as 

VA; e N, Vt > 0, j t {H k (t)e- Ct ) < H k+1 (t)e~ ct , 
thus, integrating in t, we obtain: 

VA; G N, Vt > 0, if fc (t) < / e c( ^ s) ^ +1 (s)c/s. 

Jo 

A simple recursion allows to conclude. □ 

End of the proof of prop (12.41) : Applying the previous lemma, we obtain 

Vt > 0, H (t) < — -^-y j\t - sf-'e^^'^H^ds, 
and thus, for all t G [0, T], 

i r* N 

U v - 1J- Jo i=0 

JV 



e 

< - 



(2|kll i2(/) -A)T ,t \ 



(iv-i)i y 

which gives 



i=0 



T 

sup \\fio(t, .)\\ L , m < C— x \\d" Hih^Li, (4.4) 

0<i<T ^V! ^— * 
i=0 

where the constant C depends on T, A and q. This concludes the proof of 
Prop. El □ 
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4.3 Proof of Th. 12.3b Error estimate in the general 
case 



The general idea for the proof is to regularize the initial data. Then we 
have to bound 3 terms. First the error between the exact solution and the 
truncated hierarchy for this regularized initial data. This term is controlled 
by Prop. 12.41 The next term is the difference between the solution for the 
non regularized initial data and the solution for the regularized one, both 
solutions to the exact equation (jl.4p . This is bounded using the estimates 
for (II. 4p . The final term is the difference between the solution for the non 
regularized initial data and the solution for the regularized one, but both 
solutions to the truncated hierarchy (I2.17p . We bound this term thanks to 
assumption (I2.18p . 

• Step 1: Regularization of the initial data 

We fix e > and we choose /° G H k (M. 2 ), with supp v /° C /, and such 
that (see the appendix for more details) 

a.e.vel, /°(.,d) G C°°(R), 

ll/M - f»\\ L > m < C k e k \\f(v)\\ Hkm , (4.5) 
||/ £ »||^( K ) < C k e k - m \\f°(v)\\ HHR) Vm > k. 

Let f e = f £ {t, x, v) be the unique solution to Eq. (11. 4p with f E (t = 0) = 
/°. We define the moments as usual 

fi{ s (t,x) = J ' v' ' f £ (t,x,v) dv, i G N. 

Of course, those moments still satisfy the hierarchy 

dty.{' + o x ^ +1 = 7it 4< - xtf i g n. 

We denote by \x\ the solution to the truncated hierarchy (I2.17P for the 
initial data 

H e .(t = 0, x )= fif(x) = J rffifav) dv, < i < N. (4.6) 
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Step 2: Error estimates in term of e. First note that /if — /i{ £ satisfies 
the assumptions of Prop. 12.41 On the one hand 

E n^(tf - a^ii* ^ E (n^iU 2 + iifiWi*) 

i<iV i<iV 

<En^iU 2 + eCt ( 2iV + 1 ) 1/2 ii^^ iU 2 ' 

i<iV 

by Corollary 14.31 On the other hand by applying (12. f 8p to /if and 0, 

one gets 

En^ii^<^^ 7 ii«iuiv 

i<N 

So applying Prop. 12. 4[ we have for 7' = max(7, 1/2) 



_ ^0 e |U°°([0,T],L 2 (K)) < C — ^ X ||^/ e || L 2 (K 2) 

r k-N iv 1 v 11 jrO|| 

- fc iV! 



(4.7) 



by (0311. 

We can use the stability estimate ( 12.18)) to control 

||^-/xo|U-(to^],^(M)) < CiV^||/ -/ £ °|U 2 < C fe ^ £ fe ||/°||^ iS , (4.8) 
again by (14. 5p . 

' 3 ([0,T],L 2 (R)) 



At last, we can control ||/4 — /4 £ ||l^([o,ti,l 2 (r)) according to Corollary 



11/4 ~ Mo e |U°°([0,T],L 2 (K)) - ||/4 /e |U°°([0,T],L 2 (IR)) < - / £ °I|l 2 ( 

< Cfc£ fc ||/°||^L2. 

(4.9) 

3: Choice of the parameter e 
We deduce from (14.7ft . ( 14.8ft and ( 14.9ft the complete error estimate 

/ TT2V \ 

llA*o - ^Hl-M^w) < C ^'—e k - N + N^e k j \\f \\ H , Lh 

(4.10) 
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where the numerical constant C > depends on k, T, q, and A. 

We of course choose the "best" value of e, which minimizes the error. 
When N > k, we get for 7 > 1/2 



e = e 

N\ J N->oo N 

If 7 < 1/2 then one takes instead 

' T N jy I/2-7 \ 1/N e T 



e = e 



N\ J N^oo N 



with the same asymptotic behaviour. 
In both cases 

II A*o - ^oIU°°([o,t],l 2 (r)) < Cx^xkN 1 (e*) k 

C T , qXk (4.H) 
n^oo N k -~/ ' 

which concludes the proof of Theorem 12.31 □ 



5 Proof of Th. II. lb Example of the Tcheby- 
chev points 

We now make the assumptions in Th. 11.11 and in particular that d = and 
that the satisfy (j!.16p . 

The constant A^o defined by (12. 9p satisfies 

A^ = (iV + l) 1/2 . (5.1) 

Moreover we define 



Pn[v) 



7T / 



N + l 
dizatioi 

As the Afc gives the usual method of integration we have 



In fact, the function is a normalization of the Tchebychev weight. 
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Proposition 5.1 We have, for all N > 1 and for all R £ M 2 iv+i[^] •' 



(5.2) 



We may therefore apply Th. 12.11 to this choice of and for this choice of 
Pn- 

Compute 

C N {q) = ^NTl^j\{v)\ 2 p N {v)dv^ ' <C\\q\\v. 

Similarly 

So by Th. O 

sup||^(t)|| L2 <Ce^ll^|l/°IU 2 iV- 1/2 fe A ^ ' <Ce CT ^\\f°\\ L ,, 
which is exactly f ll . 18j) . Moreover 

EW)H^^ c 'ii/ iU 2iV_1/2 Efe A ^ ' ■ 

i<N i<N \k<N / 

Of course by Cauchy-Schwarz, we have, for all < L < N, 

/ \ V2 / \ 1/2 

^ 1/2 E E a? < E E ^ 

i<N \k=L...N-L / \i<N k=L...N-L / 

I „ , V' 2 

< 



y — 

A- . i - 



\L<k<N—L 



2 



Note that 

1-|A,| 2 >%^ if ^ < iV/2, l-|A fc | 2 > (iV -^ 1)2 Xk>N/2. 
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Hence, if L < N/2, 

i<Ar \k=L...N-L / \L<fc<AT/2 / 

Therefore 

£ < g ||/°|U 2 (iv 1 ^ L V2 + iVL -i/ 2)) 

i<Af 

and choosing L = y/~N we obtain that this method satisfies the estimate 
(I2.18P with 7 = 3/4. It only remains to apply Th. 12.31 to conclude. □ 



6 Appendix 

The natural way to regularize is by convolution. However to obtain high 
order approximation, it is necessary to choose correctly the mollifier. In the 
L 2 framework though, things are quite simple by truncating in Fourier. 

Proposition 6.1 Let k a positive integer, f G H h (M. d ) and e > 0. It exists 
f £ G H°°(R d ) such that 

\\f-fe\\mRJ)<e k \\D k f\\ L2{Rd) , (6.1) 

\\D m f E \\ L2{Rd) < e k - m \\D k f\\ LHRd) , Vm G N. (6.2) 
Proof : We use Fourier's analysis. We consider f £ G L 2 (R d ) defined by 

Xco = /(0%i<im- 

First, we have 

^ ^ / r - \ 1/2 / r ~ ^ 1/2 

\\f-fe\\ L > m = i/cow < e 2k \e k \m\ 2 dt 



V|£l>l/e 

= s k \\D^f\\ LH 

which proves (I6.ip . 

On the other hand, for all m G N, 

~ / r ~ \ 1/2 / r - \ 1 / 2 

mrfeh^) = iei 2m i/(oi 2 ^ ) < ( / ^ k - 2m \e k \f{i)M 

= e k ~ m \\D*~f\\ L * { 
thus, f £ G H m (R d ) and the estimate holds. □ 
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